Neutrinos and Nucleosynthesis in Supernova * 



U. SolisQ and J. C. D'OlivcH 

Instituto de Ciencias Nucleares, 
Departamento de Fisica de Altas Energias, 
Universidad Nacional Autonoma de Mexico (ICN-UNAM). 
Apartado Postal 70-543, 04510 Mexico, D.F., Mexico. 

Luis G. Cabral-RosettH 

Departamento de Posgrado, 
Centro Interdisciplinario de Investigation y Docencia en Educacion Tecnica ( CIIDET), 
Av. Universidad 282 Pte., Col. Centro, 
A. Postal 752, C. P. 76000, 
Santiago de Queretaro, Qro., Mexico. 

The type II supernova is considered as a candidate site for the production of heavy elements. 
The nucleosynthesis occurs in an intense neutrino flux, we calculate the electron fraction in this 
environment. 
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I. NUCLEOSYNTHESIS IN SUPERNOVA 

A star lives a luminous life by burning H into succes- 
sively heavier elements. However, as the Fe group nuclei 
near mass number A = 56 are most tightly bound, no 
more nuclear binding energy can be released to power 
the star by burning Fe. Therefore, heavy elements be- 
yond Fe have to be made by process other than normal 
stellar burning. One such process is the rapid neutron 
capture process, or the r-process. 

One starts with some nuclei and lots of neutrons, 
the nuclei rapidly capture these neutrons to make very 
neutron-rich unstable progenitor nuclei. After neutron 
capture stops, the progenitor nuclei successively beta- 
decay towards stability and become the r-process nuclei 
observed in nature. This process is responsible for ap- 
proximately half the natural abundance of nuclei with 
mass number A > 100 Q. 

There is as yet no consensus for the site or sites of 
r-process nucleosynthesis. The high neutron densities 
10 20 cm 3 and temperatures of 10 9 K associated with r- 
process suggest astrophysical sites as core-collapse type 
II or lb supernovae. The most plausible environment yet 
proposed is the neutrino-heated ejecta from the nascent 
neutron star. Close to the neutron star, the temperature 
is several Mev and the atmosphere is essentially dissoci- 
ated into neutrons and protons. As the neutrinos emitted 
by the neutron star free-stream through this atmosphere, 
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some of the v e and V e are captured by the nucleons and 
their energy is deposited in the atmosphere. The atmo- 
sphere is heated and, as a result, it expands away from 
the neutron star and eventually develops into a mass out- 
flow, a neutrino- driven wind |2|. After the shock wave has 
propagated out to several hundred kilometers, condition 
behind the shock at 100 to 200 km are suitable for neu- 
trino heating. The neutrino heating blows a hot bubble 
above the proton-neutron star. 



II. NEUTRINO-NUCLEON INTERACTION 

Neutrinos and antineutrinos of all three flavors are 
emitted by the neutron star producided in a supernova. 
The individual neutrino species has approximately the 
same luminosity but very diferent average energy. As the 
neutrinos diffuse out of the neutron star, they thermally 
decouple from the neutron star matter at different radii 
due to the diference in their ability to exchange energy 
with such matter. Neutrinos species of all flavors have 
identical neutral current interactions but, due to energy 
threshold effects, the v^, v r , and their antiparticles lack 
the charged current capture reactions analogous to 



n — > p + e 
-^p + V e . 



(1) 
(2) 
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The result is that v T , v^, and their antiparticles, have 
indentical spectra and decouple at a higher density, and 
thus temperature, than the electron neutrinos and an- 
tineutrinos 0- 

We require the neutron-to-nuclei ratio R > 100 to ef- 
fect a good r-process yield for heavy nuclear species, but 
the models with conventional equations of state for nu- 
clear matter all give smaller values of R. In all mod- 
els R is determined by the net neutron-proton ratio; the 
entropy-per-baryon in the ejecta, s; and the dynamic ex- 
mx pansion time scale, td yn ■ The neutron-to-proton ratio 



2 



is n/p — Y~ l — 1, where Y e is the number of electrons 
per baryon. The r-process is only possible when Y e < 0.5 
at freeze-out from nuclear statistical equilibrium. The 
value of Y e in the region above the neutrinosphere is de- 
termined by the interactions in Eq. (1) and (2). We can 
write the rate of change of Y e with time as 

^=X 1 -X 2 Y e , (3) 

where Ai = X nl , e + A e +„ and A2 = Ai + A„„- + X ep , are 
the rates in (1) and (2), (see Ref0|). 



is the distribution function of the neutrinos or elec- 
trons and positrons. (1 — f e ) and (1 — /„) are the 
Pauli blocking factor and r is evaluated in 10 km. and 
A = ^r(gv 2 + 3gA 2 )- G F is the Fermi coupling con- 
stant, and P (Q) is a factor of survival probability (< 1) 
for neutrinos (antineutrinos) . In this work, we assume 
the energy spectrum of each neutrino species, electron 
and positron is approximated by a Maxwell-Boltzmann 
distribution, following Bernstein ||| and Enqvist Q. 

We changed to the variable y = Am/T and the solu- 
tion is function of the temperature, the integrating factor 
now becomes 



III. ELECTRON FRACTION 



The general solution to the above equation is given by 



± ( \i{t) 
dt' ^AatO, 

with the integrating factor given by 



A 2 Jo 



I(t,t) = Exp - / X 2 (t")dt 



dt , (4) 



(5) 



and A1/A2 = Y eq . The functional form for the weak- 
interaction rates are given by: 
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where B = 9.6 x lQ- u cm 2 /MeV 2 , L u = lx 10 51 erg/seg 
and Lp = 1.3 x 10 51 erg/seg are the neutrino luminos- 
ity < E v >= 3.15 T v , F a (0) = 1.8, / - E 2 /{e E / T + 1) 



i(y, y ) = Ex p 



dt 
dy" 



My"W' 



(10) 



To evaluate the integrating factor we assume the 
neutrino-driven wind model, the radius of an out- 
flowing mass element is related to time t by r = 
r exp {{t — t )/tdy n } this implies an outflow velocity pro- 
portional to the radius, v = r/td yn 7]. In the model, 
par- 3 and Tar" 1 . We set T n = T p ~ = T e = T e+ = T, 
but we assume that the neutrino and antineutrino tem- 
peratures remain nearly constant. We finally compute 
the integrals in the solution for a range of y values, we 
obtain the curve for Y e (y) and P. Numerical supernova 
neutrino transport calculations show that Tp e — 5.1 Mev, 
T Ue = 4.5 Mev and we set for Td yn — 0.3 s. Neutrino os- 
cillations add a new complication to the diagnostic of 
supernova neutrinos, the factor P will be the survival 
probability of the electron neutrinos and antineutrinos. 



IV. CONCLUSIONS 

In Figures (a) and (b) show that difference betwenn Ye 
(doted line) and Yeq (solid line) is very little. That is, the 
integral in the equation (4) not contribute significantly 
to Ye. Canculations made for Td yn between 0.01 s and 
1 s confirm the last. 
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FIG. 1: In (a) the electron fraction Y e (doted line) and Y eq (solid line) is plotted against y = (m n — m p )/T with P — Q = 
0.2, 0.4, 0.6, 0.8, 1.0. In (b) the same but assuming that no exist antineutrino oscillations, i. e. Q = 1. 
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